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A VECTOR STUDY OF LINEARIZED SUPERSONIC FLOW 
APPLICATIONS TO NONPLANAR PROBLEMS! 


By Joun C. MARTIN 


SUMMARY 


A vector study of the partial-differential equation of steady 
linearized supersonic flow is presented. General expressions, 
which relate the velocity potential in the stream to the conditions 
on the disturbing surfaces, are derived. In connection with 
these general expressions the concept of the finite part of an 
integral is discussed. 

A discussion of problems dealing with planar bodies is given 
and the conditions for the solution to be unique are investigated. 


Problems concerning nonplanar systems are investigated, and — 


methods are derived for the solution of some simple nonplanar 
bodies. The surface pressure distribution and the damping in 
roll are found for rolling tails consisting of four, six, and eight 
rectangular fins for the Mach number range where the region of 
interference between adjacent fins does not affect the fin tips. 


INTRODUCTION 


In the presentation of the theory of the flow of an idealized 
incompressible fluid, vector methods can be used to reduce 
greatly the mathematical manipulations involved. The 
study of steady linearized supersonic flow may also be aided 
by the use of vector methods. Two types of approaches, 
however, can be used. Perhaps the more obvious is to make 
use of common vector methods as was done in reference 1. 
The other vector method, which was introduced by Robinson 
in reference 2 and is used in this report, appears to be more 
suited to the study of the linearized partial-differential equa- 
tion of steady supersonic flow. This method allows a deriva- 
tion of a hyperbolic scalar potential and a hyperbolic vector 
potential along lines analogous to the derivation sometimes 
used (ref. 3, ch. VIII) in dealing with common scalar and 
vector potentials. 

The present report presents a vector derivation of many 
general results which have been found by various methods 
and are given in the published literature on the linearized 
partial-differential equation of supersonic flow and also 
presents some results which are not found in the literature. 
The general results of Hadamard (ref. 4, p. 207), Puckett 
(ref. 5), and Heaslet and Lomax (ref. 6) are found as special 
cases of a general expression for a scalar potential, and the 
results found by Robinson (ref. 2) are obtained by the use 
of a vector potential. The derivation of the scalar potential 
doubtlessly helps to clarify the concept of the finite part.of 
an integral. 


A discussion of problems dealing with planar bodies im- 


- mersed in a supersonic flow is given, and the conditions 


necessary for the solution to be unique are investigated. 
Problems dealing with nonplanar systems are also dis- 
cussed, and methods are derived for the solution of some 
simple problems dealing with nonplanar bodies. The surface 
pressure distribution, the spanwise loading, and the damping 
in roll are found for rolling tails consisting of four, six, and 
eight rectangular fins for the Mach number range where the 
region of interference between adjacent fins does not affect 
the fin tips. : i 


SYMBOLS 

A hyperbolic vector potential 

A aspect ratio of tail fin 

a positive constant 

Ci, Ca, C3 arbitrary constants 

c chord 

E, F _ arbitrary vector functions 

f scalar function defined by equation (19) 

G vector function associated with vector 
function F l 

H vorticity vector 

6/2 span of tail fin 

i,j, k unit vectors in x*-, y-, and 2-directions, 
respectively 

M Mach number 

n unit vector normal to element of area da 


ny= —if’r +jratkys 
ny = — iy! +jva' + kv 
nat = — ibn * +j" + kr” 


AC; pressure-difference coefficient 

p rate of roll 

Q function used in equation of surface of 
discontinuity 

q part of velocity vector which is made up of 

hyperbolic curl of vector potential 

qi total perturbation velocity : 

R=4V-(:—E*—PBy—n—P*"(2- 0" 

R’ small constant 

r= (2—£)*?+6t(y—n)*+- BY(2—$) 

S ` area of tail fin 

So surface of discontinuity 


1 8upersedes NACA TN 2641, “A Vector Study of Linearized Supersonic Flow Applications to Non planar Problems” by John O. Martin, 1952. 
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a a Ss, surfaces of integration 

V free-stream velocity 

Do, 01, Da volumes of integration 

1 1 

W=>z Vo—oV 5 . 

T, Y, Z Cartesian coordinates (z-axis parallel to 
free-stream direction) 

p=~+~AC—1 

T . spanwise circulation 

T. E. V (TE. 
(ef $z de=% | AC, dz) 
L. E. L. E. 

€ small positive quantity 

Ent Cartesian coordinates (axis parallel to 
free-stream direction) A 

6, p polar coordinates 

A, Y scalar functions 

A given volume 

Vi, Va, V3 direction cosines of outward normal to 
element of area da 

vy, va, Ya «direction cosines of normal (directed away 


from point (x,y,2)) to surface So 
+, n", v3* direction cosines of normal to element of 
area da used in equation (45) 


g slope of deflected area 
T area of integration 
$, $o; $1; $2, Yo scalar potentials 
C, rolling-moment coefficient per fin, 
Rolling moment per fin 
1: b 
283 
oc; 
C= » 2/2) 
V Ipo 
$ i indicates integration over closed line or 
surface 
f 
f denotes finite part of integral 
THEORY 


This report deals with the linearized partial-differential 
equation of steady supersonic flow. This equation is given 
by 

—Bdbret+bytbo2=0 (1) 


The potential is assumed to be continuous in the stream 
direction, and the potential is assumed to be always finite. 
Assuming the potential to be finite and continuous in the 
stream direction has the effect of requiring the aerodynamic 
lift and moment (calculated by use of the linearized pressure) 
of finite bodies to be finite since the linearized pressure is 
related to the derivative of the potential in the stream 
direction. The expression “linearized pressure” refers to-the 
pressure obtained by neglecting all powers of the perturbation- 
velocity components above the first. 
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VECTOR OPERATORS AND IDENTITIES 


Certain operators, which are closely associated with the 
linearized hyperbolic partial-differential equation of super- 
sonic flow (the two-dimensional wave equation), are added 
to the vector operators commonly used. The basic operators 
have been used previously in references 2 and 7. 


The gradient operator is defined by 
sO að Ò 
Y=t 3 1 g" Sz 


The analogous hyperbolic gradient operator defined by 
Robinson in reference 2 may be expressed as 


Vim ip oti tk 
The hyperbolic divergence of an arbitrary vector Æ is given by 
Vh-E 
Similarly, the hyperbolic curl of the vector E is given by 
VAXE 


The divergence of the gradient operator is sometimes denoted 
by 
D 0 


D? 
TES 


The analogous divergence of the hyperbolic gradient operator 
is denoted by | 


d a, D 
= Ve = — 8? — -L 
VhaVVi=— 6 atatoa 


The following identities are needed. Let E be a vector and 


y and A be scalar functions of r, y, and z. Then, 
Vy-VA=VAVyY (20) 
V-~E=W-E+E-Vw (2b) 
VX(VXE)=V(V-E) —V*E (2c) 
V-(VXBE)=0 (2d) 
Vhw-VA=VAA-Vy (20) 
VhYyE=yYVh-E+-E-Vhy (2f) 
VX(VRXE)=Vh(V-E) —V*hE (28) 
VAX(VXE)=V(Vh-E) —V'hE (2h) 
Vh- (VAX E)=0 (21) 


These identities can be proved by direct expansion. 
The divergence theorem may be expressed as 


| È En da =| vE dv (3) 
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where n is the normal unit vector to the element of area da. 


The vector n is expressed mathematically as 
N=, +Jretkys 


where 1, va, and v; are the direction cosines of the outward 
drawn normal to the element of area da. 

A theorem more general than the divergence theorem is 
given by (this theorem follows from the results of ref. 8, p. 
87) 


d (Ci Pi e+ Cova E,+ Cava Eda 
=| (o ta O E) de 


where the subscripts x, y, and z refer to components of the 
vector E, and Ci, Cz, and Q; are arbitrary constants. Note 
that if Ge C,=C;=1 the preceding equation reduces to 
equation (3). If ` 


G= —f* è 
C=0=1 
the preceding equation reduces to 
d (— B?», E, +» E, tado (e DE, = 32) 
or 
pen, da= | Vi-E dv (4) 
where 


n= — 18%, tjr t kr 


If the divergence theorem as expressed by equation (3) 
is applied to a volume throughout which 


V-E=0 
then the surface integral over the bounding surface is 


pen da=0 


provided that no surfaces exist inside the volume of integra- 
tion across which the normal component of E is discontinuous. 
Similarly, if equation (4) is applied to a volume throughout 
which 

Vh-E=0 


then the surface integral over the bounding surface is 
pin da=0 


provided that there are no surfaces inside the volume of 
integration across which E-n, is discontinuous. It is in- 


847 


teresting to note, however, that surfaces exist inside the 
volume of integration across which E-n can be discontinuous 
while at the same time E-n, remains continuous. It follows 
that for such a surface n and n, must satisfy the relation 

| n-n,=0 (5) 
Let Q(z,y,z)=0 be the equation of such a surface. Then, 

| Y 

VOF FO TQ: 
and 
a , 
Ny, =- VhQ 
i Y Q+ Q+ Q 


where the subscripts indicate differentiation. Substituting 
the preceding expressions for n and n, into equation (5) yields 


. 3 og 2 2 og Bes 

OQ o 
Any solution of equation (6) set equal to zero is the equation 
of a surface across which V.E may be discontinuous while 
V,-E remains continuous. The fact that the Mach cone 
from any arbitrary point satisfies equation (6) can be easily 
verified. The equation of the envelope of the Mach cones from 
an arbitrary line also satisfies equation (6) (ref. 9, p. 106). 

FINITE PART OF INTEGRALS WHICH ARISE 
IN STEADY SUPERSONIC FLOW 

In the following sections use is made of the concept of the 
finite part of an infinite integral. This concept was intro- 
duced by Hadamard (ref. 4) and has been used by a number 
of other investigators. The concept of the finite part is, 
however, sometimes confusing. This section was therefore 
included in an attempt to give a realistic picture of the 
finite-part concept and also to present the first steps of the 
derivation of the scalar and vector potentials. 

The concept of the finite part of double integrals as 
defined by Hadamard and used in this report is different 
from the concept of the finite part of double integrals as 
defined in reference 10. The essential difference between 
these two definitions lies in the manner in which the singular 
pomts along the Mach cone are treated. 

In reference 3, page 183, a vector function is used in the 
derivation of the common scalar and vector potentials. 
The analogous vector function based on equation (1) is 


1 1 


R=/(:—EP—PBly—n Be)? > 


The hyperbolic divergence of vector W with respect to 
variables £, y, and ¢ is given by 


where 


ido 9,0%, 
VA =p =p (A 
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The preceding equation indicates that the hyperbolic diver- 
gence of the vector W set equal to zero yields the partial- 
differential equation of linearized supersonic flow. A 
mathematical derivation of W can be obtained; however, 
for the purposes of this report such a derivation is not 
needed. 

The result of applying equation (4) to the vector W is 


1 1 1 
Je Vo—@V 5m da= | Vko dv (7) 


When ¢ satisfies equation (1) throughout the volume of 
integration, the right-hand side of equation (7) is zero; 


thus, 
1 1 
Je Vo— $V zm da=0 (8) 


Vihg9=0 


when 


Equation (7) is applied to a volume (denoted by vo) 
enclosed in the forward Mach cone from the point (z,y,z). 
This volume is bounded by the surface given by R=R', 
where R’ is a small constant, and an arbitrary surface S; 
enclosed in the forward Mach cone from the point (2,y,2). 
A cross section of the region of integration is shown in 
figure 1. Note that this region is analogous to the region 
that is sometimes used in calculating the potential function 
satisfying Laplace's equation (ref. 3, pp. 151-153). For 
regions such as the one shown in figure 1, equation (7) may 
be written as 


re $v t) dat Í (qv — dV: Din de 


= lz 1 vhedv (9) 


(x, ¥,2) 


FIGURB 1.—Cross section of the region of integration used in connection 
with equation (9). 
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where T represents the area of integration when R=’. 
The integral over the area T may be reduced to 


f a = (ve m2) da (10) 


where r is given by 


r= @— P+ E + eS 


Since R' is a constant, equation (10) can be written as 


y f (vent E?) da (11) 


Equation (9) can now be written as 


-4 [,(v6-m+22) dat RE Vé—87 p) 7 da 
=f J V'he dv (12) 


If ¢ is required to satisfy the linearized partial-differential 
equation of steady supersonic flow, then 


Vtho=0 
and equation (12) reduces to 


mb, (vom +42) da+ |. (5 Vp—pV 7) m da=0 (13) 


If R’ is made smaller and smaller the integrand of the integral 
over the area T'in equation (13) remains finite except-on the 
small area close to the point (x,y,2). In anticipation of tak- 
ing the limit of equation (13) as R’ approaches zero, the small 
area close to the point (z,y,2) is removed from the area T. 
The area Tis divided into two parts. One part is the area of 
T which is downstream of the surface given by 


E=t—E€ 


where eis small but larger than R’. This area is denoted by r. 
The remaining part of T (denoted by T”) is the area of T 
which is upstream of the surface 


=r— e 


A cross section of the region of integration with T divided 
into 7 and T” is shown in figure 2. Equation (13) can now be 
expressed as 


1 B4 1 | 1 1 
y | (vent ES) datiy |, (7 Vod— Y z) Ny da-+ 


1 1 
fo € Vo—oV z) -n, da=0 (14) 


where R’ is smaller than e. 

Since ¢ is continuous and therefore its values over r are 
approximately constant for small values of e, the integral 
over the area 7 can be written as 


B76 (z,y,2) f da 
val vé", dat | = | (15) 
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FIGURE 2.—Cross section of the region of integration used in connection 
with equation (14). 


When the second integral of expression (15) is integrated, 
equation (14) becomes 


yS Ven dat ENP oro (x,y,2)-+ 


L Bid f é ieee = 
1 f (Vent : ) da , (q V4—ev p) Jazi 


(16) 


If R’ is made to approach zero, equation (16) applies even 
to the limit where R’ is zero. 

The limit of equation (16) as R’ approaches zero may be 
written as 


A 2rep(1,Y,2) 
lim Fa Voy da+ R’ 219 (2,y,2)+ 


ere f G _evi ES 
Flue da+ A p V$ ¿Y p) m da i=0 (17) 


The integrands of the integrals over the areas 7 and 7” are 
always finite and it can be shown that their first derivatives 
with respect to R’ approach zero as R’ approaches zero; there- 
fore, the product of 1/R’ and these integrals either approaches 
zero in at least the order of R’ or approaches infinity as R’ 
. approaches zero. Thus it follows that the integrals over the 
arcas r and T” have no finite terms remaining after the limit 
(R'> 0) has been taken. The sum of the terms of equation 
(17) must be zero; thus the singularities resulting from the 
integrals over the areas r and T’ must cancel the singularities 
which arise from the integral over the area Sj. 

From the preceding considerations it follows that one 
method of evaluating the finite part of infinite integrals of the 
type appearing in equation (17) is to evaluate the integral 
when K’ is small but not zero and neglect the terms multi- 

321095—b5——54 
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plied by powers of 1/R’. Other infinite integrals sometimes 
arise, however, for which the finite part cannot be obtained 
by neglecting powers of 1/R”. For example, if equation (17) 
is differentiated with respect to one of the variables (x, y, or 2) 
an equation containing the velocity component is obtained. 
In some cases, when the point (x,y,z) lies on the surface S, 
the infinite terms are of the order (In R’)/R’ and of the orders 
(1/R’)*. In these cases, the finite part of the infinite integrals 
can be obtained by evaluating the integrals when R’ is small 
and neglecting the terms multiplied by powers of 1/R” and 
(In R/R. 

The process of removing the infinite parts of an integral, 
however, has been derived by Hadamard (ref. 4, book IL, 
ch. I). Hadamard used his methods of evaluating the finite 
part of integrals in finding solutions to certain hyperbolic 
equations including the linearized equation of steady super- 
sonic flow. Perhaps a fact worth noting is that the integrals 
of equation (17) are double integrals and when the methods 
given by Hadamard are used the methods given for multiple 
integrals should be used. In the past, the smgular points 
(points on the Mach cone where the derivative of 
(2— £)?— 6? (y—7)?—B?(2—f)? with respect to the variable 
of integration is zero) have caused some confusion; as 
Hadamard points out (ref. 4, p. 147), these singular points 
must be removed from the area of integration before the 
finite part is taken. Particular attention should be given to 
paragraph 92 of reference 4 since the special type of integrals 
discussed therein sometimes arises in dealing with planar 
problems. 

Robinson (ref. 2) has shown that when using Hadamard’s 
methods the order of integration may be changed without 
affecting the finite part and that it is permissible to differen- 
tiate under the integral sign of a multiple integral without 
considering the variable limits which lie along the boundary 
where the integrand is singular, provided that only the finite 
part is taken. Both Hadamard and Robinson have shown 
that in differentiating an improper integral which has an 
integrand that has a one-half power singularity along 
variable limits the variable limits may be neglected provided 
the finite part of the resulting integral is taken. 

The term “finite part” is somewhat misleading since the 
finite part of an integral can be infinite. In certain cases the 
integral is infinite even after the terms which approach 
infinity as Æ’ approaches zero have been neglected. 


SCALAR POTENTIAL 


The preceding arguments show that the finite parts of 
equation (17) can be equated to zero; thus, 


f 1 1 
| 0,4 |, (4V9—9Y5) -n, da=0 
where the symbol f before the integral denotes that only the 
finite part is to be taken. ‘The preceding equation may be 


solved for the value of the potential at the point (x,y,z); 
the result of this operation is given by 


17 1 1 
$(2,Y,2)=5— fo (g2) - Nn, da 


It should be remembered that surfaces can exist inside the 


- (18) 


— — - — -ma o aM a$ - - å; — — 
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forward Mach cone from the point (z,y,z) across which 


aVe—dV5 can be discontinuous and across which 


R 
(z Vọ—¢V z) - Na remains continuous. 


Equation (18) is an expression for the scalar potential at 
the point (x,y,2) in terms of the potential and its derivatives 
with respect to n, on the surface ©. A more general 
expression for the scalar potential than that given by 


equation (18) can be obtained. If within the volume. 


enclosed by the forward Mach cone from the point (x,y,z) 
and the surface S, equation (1) is not satisfied and V*A 


operating upon ¢ yields 
Vhe (£,9,5) =f ($m, t) 


then equation (10) becomes 


Lf (co- melt) act f, (Gro-ovi)-mäa 
= f End dv ` (20) 


Provided that /(£,9,7) is always finite, the right-hand side of 
equation (20) is finite; furthermore, the right-hand side of 
equation (20) remains finite as R’ approaches zero. If in 
equation (20) R’ is made to approach zero and only the finite 
parts of the integrals are retained, then the resulting expres- 
sion is 


J 
ó(2,9,2)=>> fa (gY$—ov 5) - roda —3> -f HEE do 


(21) 


(19) 


where v, represents the volume v when R’ is equal to zero. 
Equation (21) is equation (58) of reference 4 where f? has 
been set equal to one. Note that the volume integral in 
equation (21) has the appearance of the integral for the 
potential resulting from a volume distribution of sources 
in an incompressible flow. 

The assumption has been made that œ is continuous 
throughout the volume v, It is also assumed that no 
surfaces exist inside v, across which ò¢/òn, is discontinuous. 
If equation (21) is applied to a volume v,, which has surfaces 
across which ¢ and/or the derivative of ¢ in the direction 
of n, is discontinuous, these surfaces of discontinuity can 


be removed from the volume of integration by allowing the- 


arbitrary surface S, to envelop them (see fig. 3). For vol- 
umes of integration where the surfaces of discontinuity 
have been removed in this manner, the scalar potential 
can be written as 


aya - f LELY ap + = (hi$ 


08577 5) ds Tar Sa, 5 (BI ES n, da 
(22) 
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where Sọ denotes the surface of discontinuity, and Ag is 
the potential difference across the surface Sy. The notation 
0/dn,’ is used to denote the operator 


Ò Ò Ò 
— q? —+ ? + re 
: B n’ Of V2 On V3 Of 


For the cases where no surface of discontinuity exists 
inside the volume v; and ¢ and V¢ are zero on the surface 
S¡—S0, equation (22) reduces to 


ie,y2)——ge |, ERË ay (23) 


From equation (19) 
Vh¢ (2,Y,2) =f (2,4,2) 


Note that equation (23) is a solution of the partial-differential 
equation (24). 

For most problems in linearized supersonic flow, f(£,n,f) 
is zero and ¢ is zero upstream of the disturbing body. For 
such problems, the surface S,—S) can be taken to be located 
upstream of the disturbing body where ¢ and V¢ are zero. 
In this case, equation (22) reduces to 


L 1 Ò o 1 
y $(2,Y,2)=5— in E A ny ¿— Ad Onn’ z) da (25) 


Tf the surface So is confined to the ¿=0 plane, equation (25) 
reduces to equation (10) of reference 6. In this reference 
the boundary conditions for airfoils are discussed. 


(24) 


COMPONENTS OF VECTOR FIELD 


Let F be a vector which is finite and integrable in a given 
volume (denoted by A) and is zero outside the volume A. 
To each point in the volume associate the vector 


Gílzx,y,2)= f N Eat 15) dv (26) 


where d, denotes the part of the volume A enclosed in the 
forward Mach cone from the point (z,y,2). 
From equations (24) and (26), it follows that each com- 


ponent of G satisfies the relation 
VAG (2,y,2) = —2rFy (2,y,2) (27) 


where the subscript 7 refers to any component of the vector 
field. 


Let Yo(x,y,2) be a scalar and A(z,y,2) be a vector defined 
by the equations 


Qari =Vhi-G= f, F(6,0,)-Vh $ do (28) 
and 


; 
ITA =VX G= f Fn OXY a de (29) 


Equation (2h) indicates that 
VAX (VX G)=V(VA-G)—VLG (30) 
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Z 


(%, ¥,2) 


e 


FIGURB 3.—A cross section of the region in the forward Mach cone 
from the point (x,y,3) showing the surface 8; enveloping a surface 
of discontinuity So. 


Substituting the expressions for VX G, Vh. G, and VAG 
as given by equations (27), (28), and (29), respectively, into 
equation (30) and solving for F yields 


F(x,y,2)=—Vyole, y, 2) HVRXA (zx, y, 2) (31) 


Since F is an arbitrary vector, equation (31) indicates that 
any finite integrable vector field can be expressed in terms 
of the gradient of a scalar and the hyperbolic curl of a vector. 
Equation (31) has the appearance of the Helmholtz theorem 
(ref. 3, p. 187); however, since yo and A are found by inte- 
gration only in the forward Mach cone from the point 
(x,y,z), equation (31) hardly seems to be a statement of 
the Helmholtz theorem as is commonly given. The result 
given by equation (31) was obtained by Robinson in 
reference 2. 
HYPERBOLIC VECTOR POTENTIAL 


Equation (31) indicates that the perturbation velocity 
vector can be divided into two parts. One part is the 
gradient of a scalar function, and the other is the hyperbolic 
curl of a vector function. The vector function is analogous 
to the common vector potential (ref. 3, pp. 104 and 188); 
therefore, the vector function is referred to as the hyperbolic 
vector potential. Thus, if q' is the total perturbation 
velocity vector, then 


q’=VetVAXA 
where ¢ is the scalar potential and A is the hyperbolic vector 
potential. The part of the velocity vector which is made up 


of the hyperbolic curl of the vector potential is denoted by q. 
By direct expansion it can be shown that 


Vivg’=Vh-Vo+Vh- (VAX A)=VhO+Vh-(VkXA)=0 (88) 


(32) 
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Equation (33) indicates that the hyperbolic divergence of 
the perturbation velocity vector is zero. 
The vorticity vector is given by 


H=VXq’' (34) 


Therefore, from equation (32), 


H=VX(vhXA) 
or 
H=vh(v-A)—VhA 


From equations (2d) and (29), the divergence of the hyper- 
bolic vector potential is zero; thus, 


H=—WhA (35) - 


Each component of equation (35) is a partial-differential 
equation of the form of equation (24); thus, from equation 
(23) each component of equation (35) has a solution given by 

— a. Hen, 9) 

Ade =p, | HEED do (36) 
where the subscript 2 refers to any component of the vector 
H. Since each component of A is given by equation (36), 
then 

= 1 A(E,n, 6) E 
A(z,y, =z, R dv (37) 


The velocity vector resulting from the hyperbolic vector 
potential is therefore given by 


g=vi xa = | vaxi dy (38) 
91 
or l 

27 

a = ENE 7 y 
ya a E n)H. a DH, dy (398) 
F 1 
2 f i —[y — 


where the subscripts refer to the components of the vector 
H. The results given by equations (39) were obtained by 
Robinson in reference 2. 


VORTEX SHEETS 


If the vorticity is confined to a surface S3, equation (37) 
becomes 


1 H(£,m,4) da 


A(z,y,2)=5- & R 


(40) 
Equation (40) is an expression for the hyperbolic vector 
potential resulting from a surface of vorticity. Note that 
if the vorticity is zero except on the surface Sz, then equation 
(35) reduces to 
VhA=0 


By removing the surface S} from the volume of integration 
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each component of A can be expressed as (from eq. (25)) 


a oft 
A2,y,z)= on E R A Sn 


where the subscript 1 refers to any component of the vector 


A. Since each component of A is given by equs tión (41), 
then 


1 AoA 
Ao (ps 


Note that if AA is zero, equation (42) reduces to 


Alz,y,2) = EA ga da 


eh ong (43) 


By comparing equations (40) and (43) it follows that, on the 
surface S3, 

E AO 

H= Asa? (44) 

Equation (44) indicates that across a surface of vorticity 

_ the derivative of the hyperbolic vector potential in the di- 
rection of n,’ is discontinuous. Thus, a lifting surface can 
be represented by a continuous hyperbolic vector potential, 
while it can be shown that a thickness effect can be repre- 
sented by a discontinuous hyperbolic vector potential. 
Note the contrast with the scalar potential, which uses a 
continuous potential to represent a thickness effect and a 
discontinuous potential to represent a lifting surface. 


FURTHER DEVELOPMENT OF SCALAR POTENTIAL 


The scalar potential can be expressed in forms other than 
those already presented. Equation (8) is applied to the 
region bounded by the arbitrary surface S,, the forward 
Mach cone from the point (x,y,z), and a second arbitrary 
surface S3 enclosed in the forward Mach cone from the point 
(x,y,2) and upstream of the surface S,. A cross section of 
such a region is shown in figure 4. The result of applying 
equation (8) to this region is 


¡ME pl oe 
Zar on,* 
(45) 


provided that ¢’ is a solution of equation (1). Note that the 
scalar potential as given by equation (18) is independent of 
¢’ so that $” is arbitrary so long as it satisfies equation (1) 
throughout the proper volume. . 

If for a finite distance upstream ¢’ is zero and remains zero 
for greater distances upstream, the surface S, may be chosen 
in this region so that the integral over Ss i in equation (45) is 


zero; thus, 
—4 f 1 òg 
2r S1 R ón 
Equations (18) and (46) can be combined to yield 


Op _ dp , 
$(x,y,2)= NEO On, Dm" —(ġ— Ns 


1 09 


O Ns 
PA 52s aa fG R dns —¢ 54 p) 


NS 


{do 


(47) 


—AÁ ><> On,’ 7 T da (41) : 


AA Sr 5) da (42) | 
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R=0 


(x, ¥, 2) 


FricGuRE 4.—Cross section of the region of integration used in connection 
with the potential function 9”. 


The only restrictions placed on ¢’ at this point are that it 
satisfy equation (1) and be zero at a finite distance upstream. 
dp one i 
On, On, 

is zero; therefore, in these cases, ¢ can be road as 


let òp 2, 
s R On, 


1 / o 1 
$(2,Y,2)= o if (p—9”) om, R da 


In many cases ¢’ may be chosen so tha 


$(x,y,2)= (482) 


or 
(48b) 


Equations (48) are quite useful; however, remember that 
they apply only when ¢’ can be chosen so that ¢’ does not 
violate any of its restrictions. 

Note that equations (48) can be applied to problems whore 
either $ or 0¢/On, is given on the surface S,. The application 
of these equations to most nonplanar problems of either type, 
however, lead to quite unwieldy integral equations. 


APPLICATIONS 
_ PLANAR PROBLEMS 


Many problems in linearized supersonic flow deal with the 
surface of discontinuity confined to a plane surface parallel 
to the z-axis. In this section a general discussion of this 
type of problem is given. The coordinates are located so 
that the surface of discontinuity is in the ¿=0 plane. 

The scalar potential at an arbitrary point (x,y,z) above the 
ae plane is (from eq. (47)) 


seno, [te h+0—0) EE] at dr 
(49) 


In this case, the surface S, is the ¿=0 plane. 
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If $' (x,4,0) is chosen equal to ¿(x,y,0) the potential becomes 
symmetric with respect to the ¿=0 plane. Thus, for z=4 


$ (1,y,0) = p' (2,Y, —a) 
and 


Ps (2,y,a) =— ¢,' (£Y, —a) 


For this case, equation (49) reduces to 


a | SELO de dy (50) 


Equation (50) was given by Puckett in reference 5. 

If $'(x,y;0) is chosen equal to —¢(z,y,0) the potential 
becomes antisymmetric with respect to the ¿=0 plane. In 
this case, 

$(2,y,a)=—¢’ (z,y,—@) 
and 
bz (2,y,a) = by (x,y, —a) 


Thus, equation (49) is reduced to 


ln f 
a f SENO de da (51) 


Note that for surface S, not confined to a plane parallel to the 
z-axis, 2 choice of q'(x,y,2) at the surface S; to equal ¿(x,y,z) 


at S, does not causo 222,2) at the surface S, to equal 
h $ 
/ 
Epa) at Si. Similarly, choosing ¢’ (x,y,z) at the surface 
h 


/ 
S, equal to —¢(z,y,z) at Sı does not cause See a) at Sı to 
A 


equal ce) at Si. 

eee the discontinuities are restricted to the ¿=0 
plane, the scalar potential can also be expressed as follows 
(from eq. (18)): 


rl, o FF [ae dn 
(62) 


for positive z. A comparison of equations (50) and (51) with 
equation (52) shows that the two terms of the integrand of 
equation (52) contribute equal amounts to the potential at 
any point (x,y,2). 

Since the terms of the integrand of equation (52) contribute 
equal amounts to the potential at the point (z,y,z) as z 
approaches zero, equation (52) must reduce to 


¢(2,y,0)—=— zla Lib nd) de dn pE 

The preceding equation can also be obtained by examining 
the limit of equation (52) as z approaches zero. If this 
procedure is done the entire contribution of the second term 
of the integrand of equation (52) is found to come from the 
‘ point at the apex of the hyperbola formed by the intersection 
of the Mach cone from the point (x,y,2) and the ¿=0 plane. 
Note that if the integration is performed first with respect to 
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7 then, when the aon of Hadamard are used, the point at 
the apex of the hyperbola is a singular point and must be 
removed from the area of integration by a process such as is 
given in reference 4, page 147. 

If $, is prescribed over the ¿=0 plane, then the potential 
is given uniquely by equation (50). Similarly, if the potential 
is prescribed over the ¿=0 plane, then the derivative of ¢ with 
respect to z is determined over the z=0 plane. This result 
follows from equation (51) since prescribing ¢ over the 2=0 
plane determines the potentials in the space above the 2=0 
plane; therefore, it also determines ¢, in the space above the 
z=0 plane and the limit of ¢, as z approaches zero from the 
positive direction. 

The question that arises is whether ¿(x,y,2) is uniquely 
determined in the space above the 2=0 plane if ¢ is prescribed 
over certain areas of the 2=0 plane and ¢, is prescribed over 
the remaining areas.. If the assumption is made that ¢ is 
not determined uniquely, then at least two potential functions 
satisfy the condition that either ¢ or $, is prescribed in all 
regions on the z=0 plane and that ¢ is identically zero 
upstream of a given point. Let $, and ¢ denote two poten- 
tial functions which satisfy the same boundary conditions, 
and let ¢) denote the potential function formed by taking the 
difference between ¢; and $3. Mathematically, the potential 
function de is given by 


Go (1,4, 2) =¢ (x,y, 2) — do (zY, 2) (5 3) 


Since $, and ds have the same values in certain regions in the 
z=0 plane then d is zero in these regions. Similarly, since 
061/02 and O¢/dz have the same values in the remaining 
regions of the z=0 plane, then 0g,/0z is zero in these remain- 
ing regions. The potential function ¢ọ has the boundary 
conditions that either ¿y or 0¢,/0z is zero in all regions of the 
2=0 plane and that ¿y is identically zero upstream of a given 
point. 

Consider the case a all the boundaries between the 
regions are supersonic. (The slopes of the boundaries are 
such that the component of the free stream perpendicular 
to the boundary is always supersonic.), The potential 
function ¢) can be evaluated by use of equation (50) or (51) 
for points in areas which are far enough upstream to be 
affected only by a region where ¢ or ¢, is prescribed. Forall 
points in these areas ¿y is zero as indicated by equation (50) 
or (51). It follows from equations (50) and (51) that ¢ is 
also zero inside the volume above the z=0 plane, which is 
affected by these areas alone. Thus, the volume where q 


1s identically zero has been moved downstream. ‘The same 


argument can be repeated until the complete z=0 plane has 
been covered. The preceding arguments cannot be applied 
to cases where the regions have subsonic boundaries; however, 
if it is permissible to distort the boundaries within a strip of 
infinitesimal width these subsonic edges can be converted 
into supersonic edges by replacing every element of the 
subsonic boundaries by a broken line made up of supersonic 
segments. Such a procedure is illustrated in figure 5. If 
the assumption is made that the subsonic boundaries may be 
distorted an infinitesimal amount, then q, is zero over the 
z=0 plane and also in the space above the z=0 plane. 
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(a) 
(a) Subsonic element within (b) Broken line replacing 
infinitesimal strip. line element. 


FIGURE 5.—An illustration of a method of replacing a subsonic line 
element with a broken line made up of supersonic segments. 


Equation (53) now reduces to . 
diz, 12) = $2(2,Y,2Z) 


Since $, and q. are any two potential functions with the same 
values in certain regions of the z=0 plane and with the same 
partial derivatives with respect to z in the remaining regions, 
proof has been given that only one potential function exists 
for which the potential is prescribed over certain areas in the 
z=0 plane and the partial derivative with respect to z is 
prescribed over the remaining areas. 

The boundary conditions for a zero-thickness lifting airfoil 
with a given local angle-of-attack distribution are not of the 
type discussed in the preceding paragraph. The conditions 
prescribed in the 2=0 plane for this type of problem are: 
The potential function ¢ is identically zero upstream of the 
airfoil; ¿(x,y,0) is zero except on the plan form or in the 
wake; the partial derivative of the potential with respect to 
2, ¢:, is given on the plan form; and ¢,(2,y,0) is zero in the 
wake. The preceding boundary conditions do not specify 
that ¢ or ¢, be prescribed in all regions on the 2=0 plane 
since not ¢ but ¢, is given in the wake. For airfoils which 
have trailing edges which are always supersonic, the require- 
ment that ¢ be continuous in the stream direction necessi- 
tates the potential in the wake to have the value of the 
potential at the trailing edge of the airfoil. In this case, the 
potential function is uniquely determined. For airfoils 
which have subsonic trailing edges the Kutta-Joukowski 
condition is generally applied to the trailing edges to deter- 
mine ¢ uniquely. If the assumption is made that the trail- 
ing edge can be distorted within a strip of infinitesimal width, 
then the requirement that ¢ be continuous in the stream 
direction can be used to determine ¢ uniquely. If the as- 
sumption is made that the subsonic trailing edge is distorted _ 
within the infinitesimal strip so that each segment of each 
line element of the trailing edge is always supersonic (see 
fig. 6), then ¢ is determined uniquely. It is well-known 
that for airfoils with subsonic trailing edges there are an 
infinite number of solutions which satisfy the boundary 
conditions as stated at the beginning of this paragraph. The 
preceding arguments however prove that there is but one 
solution for an airfoil which has had its subsonic edge re- 
placed by broken lines which are always supersonic. Note 
that it has not been proved that the solution obtained by 
distorting the subsonic trailing edges corresponds to the 
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Boundaries ot. 
infinitesimal strip -s 


Subsonic 
4 trailing edge 


-Distorted 
trailing edge 


(c) 
(a) Wing with subsonic edge. 


(b) Line element at a break 
in a subsonic trailing edge. 
(c) Possible distortion of a subsonic trailing edgé at a break in 
that edge. 
Figure 6.—Methods of distorting a subsonic trailing edge to determine 
the potential function uniquely. 


solution satisfying the Kutta-Joukowski condition, nor has 
it been proved that the solution of the distorted trailing edge 
is independent of the manner of distortion. 


NONPLANAR PROBLEMS 


The scalar potential resulting from the disturbances caused 
by a nonplanar body can be found from equation (18) 
provided that both ¢ and 0¢/On, are known on some surface 
Sı. Unfortunately, ¢ and 0¢/dn, are not generally known 
on a surface which fills the requirement of the surface Si; 
therefore, equation (18) appears to have little value in the 
calculation of the potential functions for nonplanar systems 
in general. Certain properties of equation (18) are, however, 
worth investigating. 

The problem of evaluating the potential on the upper sur- 
faces of a long rectangular body is discussed. The assump- 
tion is made that the body extends upstream to infinity and 
that the sides are parallel to the free-stream direction except 
for small local variations which cause small disturbances in 
the stream. Figure 7 (a) shows the forward Mach cone from 
a point on the upper surface of such a rectangular body. 
This figure also shows that there is a certain part of the sur- 
face of the rectangular body in the forward Mach cone from 
the point (x,y,z) that cannot possibly affect the potential at 
the point (x,y,z). If the surface S, in equation (18) is 
taken to be the surface of the rectangular body, then equation 
(18) indicates that the values of ¿ and 0g4/0n, in the region 
which cannot possibly affect the potential at the point 
(x,y,z) should be used in evaluating the potential at the 
point (x,y,z). The only possible explanation of this con- 
sideration is that the integral of ¢ and 0¢/On, caused by the 
disturbances in the “blind spot? add to zero. This con- 
sideration can be shown mathematically as follows. Let Yo 
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denote the potential function resulting from the disturbances 
inside the blind spot. From equation (46), it follows that 


I Ovo o 1 
¡ME s,\R dn," Yon B aa ae (52) 
Equation (54) indicates that the potential at the point 
(x,y,2) can be evaluated by applying equation (18) to the 
surface of the rectangular body regardless of the blind spots. 
The same argument holds for other bodies with blind spots. 
The preceding arguments can be clarified by a simple 
illustration of the effect of a blind spot. Consider an infinite 
rectangular body such as shown in figure 7 where the only 
disturbances are caused by a small deflected area with a 
constant slope o with respect to the free-stream direction 
located on the lower surface of the body. The leading edge 
of the deflected area is chosen perpendicular to the free- 
stream direction so that the potential in the region not 
affected by the vertical sides is of a two-dimensional nature. 
Figure 7 (b) illustrates such a disturbing surface. _ 
The disturbance potential in the two-dimensional region is 
given by 


$ (2,2 EE (2—22)] 


where the lower surface of the body lies in the z=z, plane 
and the leading edge of the deflected area i8 in the z=x* 
plane. 


Forward Mach cone 


-(4,4,2) 


Intersection of the forward -~ 
Mach cone with the 
rectongular body 


(9) 


(a) The forward Mach cone from a point on the upper surface of a 
rectangular body. 


Fiours 7.—Rectangular body parallel to free-stream direction. 
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(b) ae view of rectangular ee with deflected area on the lower 
surfa 


„~ Detlected area 


FIGURE 7.—Continued. 
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_ 


-r ‘er cs a A e 
+ 


_- Forward Mach cone 


3 r(x,y,2) 


N-- Deflected area 


Lana Intersection of the forward Mach 
cone with the side of the body © 


| (c) 3 
$) Side view of dE body and the forward Mach cone from the 
point (x,y,z). 
Figure 7.—Concluded. 


A point on the upper surface which has only the two- 
dimensional flow in its forward Mach cone is illustrated in 
figure 7 (c). The disturbance potential for this point is 
(from eq. (18)) given by 


f 2 
$(2,y,2)= i Ile Y y poe) aed 


Upon performing the indicated integrations the preceding 
expression becomes 


__oV [ex—x%2—B (2—29)] , eV [x—22—8 (2—22)] 
$(2,y,2)= E E Y ee 


which reduces to 
$ (2,y,2)= 0 


This result is a demonstration that the disturbances in blind 
spots do not contribute to the potential. 

The scalar potential resulting from the disturbance pro- 
duced by a nonplanar body can also be obtained by use of 
equations (47) and (48) provided that the necessary values 
of $, 0¢/On,, $”, and O¢’/dn, are known. In dealing with 


planar bodies ¢’ could be chosen so that 


$—9'=0 


and thus equation (47) is reduced to equation (50). 
larly, ¢’ could also be chosen so that 


Od Og” 
on, on,* 


Simi- 


and thus equation (47) is reduced to equation (51). Un- 
fortunately, for nonplanar bodies, choosing ¢’ equal to $ 
does not make 0¢’/On,* known as was the case for planar 
problems and, similarly, choosing 0¢’/dn,* equal to —d¢/dn, 
does not make ¢’ known. Certain problems exist in which 
¢’ can be chosen so that ¢ can be written as a simple integral. 


INTERSECTING PLANES 


Many problems concerning nonplanar bodies deal with 
disturbances produced by two intersecting planes parallel 
to the free-stream direction. In this section, methods of 


ee ee i Ce 
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solutions for two planes intersecting at various angles are 
'given. The component of velocity normal to the surface is 
assumed to be known. 

Perhaps the simplest case of two intersecting planes occurs 
when the planes intersect at right angles. It is desired to 
find the potential in space resulting from the disturbances 
produced by the two intersecting planes. This type of 
problem could represent an isolated cruciform tail with 


supersonic leading edges undergoing various motions. ` 


Problems of tbis type have been solved in references 11 and 12. 
The axes are chosen so that y=0 and 2=0 are the disturbing 
planes (see fig. 8). When y and z are positive, equation 


(18) becomes 

__ 1f 1266040, By 0D, _ 
py 2)= e, R On da 2r Js, IP da 

A 1 O¢(é,7,07) _ Biz! o(é,7,0*) 

Qa J Or dd 2r fa R da (55) 


The surface S; has been taken to be the disturbing surface; * 


thus, S; is the part of the y=0 plane (z positive) bounded 
by the 2=0 line in the y=0 plane and the forward Mach 
cone from the point (z, y, 2). Similarly, S, is the part of the 
2=0 plane (y positive) bounded by the y=0 line in the z=0 
plane and the trace of the forward Mach cone from the point 
(x, y, 2) (see fig. 9). 

The assumption is made that ISERO) aro ) 
known and that (¢,0*,¢) and $(&,7,0+) are unknown. The 
integrals containing ¢(é,0t,¢) and $(£,7,0+) can be eliminated by 
several applications of equation (46). Equation (46) is 
applied to the volume on the left-hand side of the y=0 plane 


oe sd 


y 


Ousturbing 
surface 
(y =0 plane) 


Disturbmg surface 


(2=0 pone) 
` A 


x 
FIGURE 8.—Two disturbing surfaces intersecting at right angles. 


RBPORT 1143— NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


enclosed by the forward Mach cone from the point (2,y,2), 
the y=0 plane, the z=0 plane, and an arbitrary surface 
upstream of the disturbance (see fig. 10). The result of 
applying equation (46) to this volume is 


Lf FG y Be" (6 0%.8) om 
3 


aJs R. 0 2r 
1 f 109'(507,7) cu oe 5,5) 
EN O A TS da a da=0 (56) 


The surface S, has been taken to bo. the y=0 plane (z posi- ` 
tive) and the 2=0 plane (y negative); thus, Ss is the part of 
the z=0 plane in the forward Mach cone from the point 
(x,y,2) (see fig. 10). Adding equations (55) and (56) yields 


Sd 
$(£,0+,7)- 09 se 7) 
p(x,y,2)=— 2 =F On ott da— 


201 
Se E- [o(é,0t+,)—9'(E,0-, Dida — 


09(£,7,0+) 
if. of 
Re Rk 
a +) 99'(€,7,0+) 
pz OE 1 dt 
mlae "ala R > 
per o'(E,0F, 8) 
a le ES (57) 


PR ee r- + — r 


(y = O plana) PS i A 


Tenna. Mach cone 
from the 
point (x,y, 2) 


FIGURE 9.—Regions of integration for equation (55). 
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The potential function ¢’ (é,7,¢) is chosen so that 
p (E, —4a,f) = $(£,a,¢) 


where a is positive. In this case, 
a¢'(z,0-,t)__ d4¢,04,5) 
On On 


and equation (57) reduces to 


` — l AO 1 1 09(£,7,0*) _ 
tol On da— O aR Òt da 
ee da— 1 109 (E,7,0*) da— 
a E O oe Jog R or 
/ + 
E ¿Egon ) da (58) 


Since ¢’(§,7,¢) is related to ¿(£,7,7) the only unknown in the 
preceding equation is ¿(£,y,04). The region of integration 
S; becomes the part of the 2=0 plane (z negative) in the 
forward Mach cone from the point (x,y,z) obtained by re- 
flecting the disturbing surface in the 2=0 plane (z Pte 
through the 2=0 plane (see fig. 10). 

The problem being considered is one in which the sauna 
derivatives of the potential function are known on two planes 
parallel to the z-axis and are intersecting at right angles. 
The point (z,y,z) has been restricted to positive values of 
y and z. For the present, consider the problem of finding 
the potential above the z=0 plane when the derivative of 


pe 


Reflection of the 
disturbing surface 
In the 2 =Q planes 


E ( V2) 
\. through the y = O planes 


point 


FraureE 10.—Regions of integration for equation (56). 
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the potential function with respect to z is known on the z=0 
plane, the derivative of the potential function with respect 
to y is discontinuous across the y=0 plane, and the potential 
is zero a finite distance upstream of a given point. From 
equation (22), the potential for this problem is given by 


genz) ==; |, q 4 SEED a — 


1 E d¢(E,n,04) da 


pz! $(£,9,0+) 
27 Js R Of Ie RS da (59) 


where S, is the part of the y=0 plane above the 2=0 line 
inside the forward Mach cone from the point (x,y,z), and S; 
is the part of the 2=0 plane inside the forward Mach cone 
from the point (z,y,z). For positive values of y and z, 
equation (59) reduces to equation (58) if 4’ (£,9,0+) is assumed 
to be the true value of ¿(£,9,0*) when y is negative. Since 
in the original problem ¿(£,n,7) was not defined when y was 
negative, nothing is violated if it is now defined as being 
$'(£,n,5) in the region where y is negative and ¢£ is positive. 
The problem in which the normal derivative of the potential 
function is known on two planes parallel to the z-axis and 
intersecting at right angles has, therefore, been changed to 
the problem in which the derivative of the potential has a 
known discontinuity across the y=0 plane (z positive) and . 
the normal derivative of the potential is known on the 2=0 
plane. Note that the potential function still remains unde- 
fined below the 2=0 plane. 

Since the potential function is undefined below the 2=0 
plane, it can now be defined so that the resulting potential 
function is symmetric with respect to the z=0 plane. De- 
fining the potential below the 2=0 plane so that 


$(E,n,a) =¢ En, —a) 


yields the desired symmetry. The result of applying equa- 
tion (46) to the region below the z=0 plane inside the forward 
Mach cone from the point (x,y,2) is 


1 1 O¢(é, UE 0*) _ bz a l, VE 0*) 
27 Js, R Of aa 2r Js W a 
1 1 o¢(é, un) 0*) aS PLE, Y) el, 7,0?) ' 
ajs B or 4 s gp t 


1 1 dolt, OF, ‘oa 
1 f oh dp da=0 (60) 
where S, is the part of the y=0 plane (z negative) inside the 
forward Mach cone from the point (x,y,2). In applying 
equation (46) to the region below the z=0 plane, the surface 
Sı had to be folded over the part of the y=0 plane (z nega- 
tive) across which 0¢/07 was discontinuous in order to be 
able to apply equation (46) to this region. Equations (58) 
and (60) can be combined to yield 


_ If 10¢G,0*,8 1f 104(£, 9, 0?) 
plz, Y, 2)= Aa On da— a R dE da — 


Lf 1266 2,09 q, 1 f 1246, 04,0) y, 
S; R TJ) So 


of 


a A nr — 
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Equation (61) contains only integrals of known expressions, 
and it is, therefore, the solution to the problem of two planes 
intersecting at right angles parallel: to the z-axis where the 
‘normal derivative of the potential function is known on both 
planes. Figure 11 shows the cross section of the distribution 
of velocity normal to the surfaces for a problem as repre- 
sented by equation (55) and its solution as given by equation 


+ 
(61). Note that a (n negative) in S; is the reflection 


+ 
of oein in S, (y positive) across the y=0 plane and 


+ + 
that O in Ss (7 negative) is the reflection of ES ib) 


in S3 (y positive) across the ¿=0 plane. This condition 
suggests that the result given by equation (61) could also be 
obtained by utilizing the concept of reflecting surfaces. 
The mathematical derivation required for finding solutions 
to problems consisting of two planes parallel to the z-axis 
intersecting at various angles can be reduced by making use 
of the concept of the reflecting surfaces. For this reason, the 
result given by equation (61) is obtained by use of reflecting 
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i (a) Original problem (equation (55)). 
_ (b) Solution to problem as given by equation (61). 
FIGURE 11.—A cross section of the distribution of the velocity com- 
ponent normal to the 2=0 and the y=0 planes represented by 
equations (55) and (61). 
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surfaces. The potential function can be separated into two 
parts, ¢, and ¿,, satisfying the following boundary conditions 
on the disturbing surfaces: 


0g1(£,91,04) _ 9g(£,9,0+) 
òt Òt = 
Ogi(E,0*,5)__ 9 
97 
O¢a(E,7,0*) o 
of 
dealt, 0+, t) _ D6(E,0+,1) 
On On 
A cross section of these boundary conditions is shown in 
figure 12. Only the potential function ¿, is treated in detail 
since the boundary conditions for $, and œ: are of the same 
type. The normal derivative of q, is zero on the y=0 


plane; thus, the y=0 plane can be considered as a reflecting 
plane. The potential function ¢, is, therefore, the potential 


+ 
resulting from a distribution of Pda which is symmetric 


+ 
with respect to the y=0 line and has the value of PE ) 


_ when y is positive. Figure 13 illustrates such a distribution. 


The problem of finding ¢, has been reduced to a planar 
problem which can be solved by use of equation (50). 
Equation (50) was obtained by defining the potential 
below the 2=0 plane so that the total potential function was 
symmetric with respect to the z=0 plane. This result 
caused the derivative of the potential function with respect 
to z to be antisymmetric with respect to the 2=0 plane. 
Figure 14 illustrates the distribution of the normal derivative 
of the function across the z=0 plane. The problem of 
evaluating the potential function ¢, has been reduced to a 
planar problem. Similarly, the problem of evaluating the 
potential function ¢ can be reduced to a planar problem. 
Figure 15 illustrates such a procedure. The original poten- 
tial function is the sum of ¢, and $3. Equation (61) follows 
from the preceding results for $, and ¢. The addition of 


~ $, and q» is illustrated by figure 16. 


The concept of reflecting surfaces is now utilized to find the 
potential resulting from two disturbing surfaces parallel to 
the x-axis and intersecting at an angle of 45°. The axis is 
chosen so that the z-axis lies along the intersection of the 
disturbing surfaces and ‘one of the disturbing surfaces lies 
in the z=0 plane (see fig. 17). The potential function ¢ is 
divided into two parts, ¢, and $3. The boundary conditions 
on ¢; and ¢; are similar to the corresponding potential func- 
tions used for the disturbing surfaces intersecting at 90°. 
Figure 18 illustrates the boundary conditions for q, and q. 
The surfaces on which the normal derivative of ¢, is zero 
can be considered as a reflecting surface. This consideration 
leads to the same distribution of the normal derivative of 
$, on the y=0 plane as is given on the ¿=0 plane. Figure 
19 illustrates such a distribution. The problem of finding : 
$, for two disturbing surfaces intersecting at 45° has been 
reduced to a problem of two surfaces intersecting at 90°. 
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FIGURE 12.—A cross section of the distribution of velocity normal to the disturbing surfaces for the potential functions ¢, $1, and q». 
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FIGURE 14,—The normal derivative of the sienta function across the z=0 plane obtained by applying equation (50) to a planar problem. 


The solution of ¢, can be obtained from equation (61). 
Figure 20 shows the surfaces across which the normal deriva- 
tive of ¢, is discontinuous. Since $, and ¢; have the same 
type of boundary conditions, then ¿3 has a solution as illu- 
trated in figure 21. The original potential function ¢ is the 
sum of ¢, and ¢; therefore, $ can be found by considering 
surfaces of discontinuity as illustrated in figure 22. The 
potential function ¢ can be evaluated by use of equation 
(22), because no surfaces across which ¢ is discontinuous 
exist and the values of A(0¢/On,) are known across all sur- 
faces of discontinuity. 

Another simple case of two disturbing surfaces parallel to 
the z-axis occurs when the surfaces intersect at an angle of 
60°. The potential function ¢ is divided into two parts, 
$, and $3. The boundary conditions on ¢; and ¢, are similar 
to the corresponding functions used previously. Figure 23 
shows a cross section of these boundary conditions. By use 
of a reflecting surface, the function ¢, can be represented by 
the boundary conditions as shown in figure 24. The func- 
tion $, is undefined for 240° of the total angle around the 
z-axis. The function q, is defined as shown in figure 25. 


Since no surfaces exist across which q, is discontinuous, the 
function ¢, can be evaluated by using equation (22). Simi- 
larly, ¢2 can be defined ds shown in figure 26. The sum of 
$, and ¢ is illustrated by figure 27. The potential function 
$ can be found by using equation (22). 

In the preceding paragraphs, methods have been found 
for determining the potential resulting from two plane dis- 
turbing surfaces parallel to the stream direction intersecting 
at certain angles. The same method can be used to find 
methods for determining the potential resulting from two 
plane disturbing surfaces intersecting at various other angles, 


ROLLING TAILS WITH MULTIPLE RECTANGULAR FINS 


The methods derived in the preceding section are used to 
find the surface velocity potential, the pressure distribution, 
and the damping in roll of rolling tails consisting of four, 
six, and eight rectangular fins. For comparison, these same 
quantities are also presented for the planar tail configura- 
tions consisting of one and two rectangular fins. An illus- 
tration of the tails treated is shown in figure 28. The. 
analysis is limited to tail configurations having surfaces of 


_FIGURE 15.—An illustration of the reduction of ¢, to a planar problem. 
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FIGURE 16.—An illustration of the addition of $, and q. 


vanishingly small thickness’ and of zero camber. The 
investigation is also limited to the range of Mach numbers 
for which the region of interference between the adjacent 
fins does not affect the fin tips (see fig. 29). 


TAIL CONSISTING OF ONE FIN 


The pressure distribution and the velocity potential on 
the surface of rolling tails made up of one and two rectangular 
fins can be obtained from the results of reference 13. The 
pressure and potential for the tail consisting of only one fin 
can be found by transforming the axis of roll of the tail 
‘consisting of two rectangular fins. 

The tail consisting of one fin is divided into regions as 
shown in figure 29 (a). The velocity potential on the sur- 


face facing the negative y-direction is given for the various 
regions as follows: 


For region I, 


plz, ¿=P |z sin”! y iA 
RE e 


PEZ 


olz, z) =- 


For region II, 
(62b) 
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For region II, 


pz, ¿=P A sin”? y D. (F+ 2) y Z G- 2) | (62c) 


For region IV (note that the potential in region IV is the 
potential in region I plus the potential in region UT minus 
the potential in region Il), 


2p | TZ in-i Pe ain- 


EEE | 
E 


(62d) 


The pressure-difference coefficient is given for the various 
regions as follows: 


p(z, 2)=— 


FIGURB 17.—Position of coordinate axes for disturbing surfaces 
intersecting at 45°. 
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FIGURE 18.—A oross section of the velocity distribution on the disturbing surfaces for the functions $, and ¢, for the 
disturbing surfaces intersecting at 45°. 


FIGURE 19.—Reflection of the normal derivative of ¢, on the y=0 plane. 


— 


862 REPORT 1143—NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


= 


FIGURB 20.—A cross section of the surfaces of discontinuity which can be used to evaluate the potential function ¢, for 
two surfaces intersecting at 45°. 


For region I, 
b 
(5—2 
slo EE 
(63a) 
For region LU, 
pz 
AC» =37 (63b) 
For region III, 
ACE E sin [Eta] 2 (5—2) | 63o 
For region IV, l 
b 
pce 
A0, =y . in VE b sin E 
x b b IN 2 (63d) 
va (5-=)YG--) ame 


TAIL CONSISTING OF TWO FINS 


The tail consisting of two fins (fig. 29 (b)) has the same 
potential and pressure distribution as a rectangular rolling 
wing and can, therefore, be obtained from the results of 
reference 13. For each tail consisting of two fins divided 
into regions as shown in figure 29 (b) the pressure and poten- 
tial in regions I and IT are the same as the pressure and po- 
tentials in the corresponding regions for tails consisting of 
one fin. 


TAIL CONSISTING OF FOUR FINS 


Each fin of the tail consisting of four fins is divided into 
regions as shown in figure 29 (c). The pressure and poten- 
tials in regions I and II are the same as the pressure and 
potentials in the corresponding regions for tails consisting of 
either one or two fins. The regions DI and IV are affected 
by the interaction between adjacent fins. The potential in 
region IV is made up of a combination of the potentials of 
regions I, II, and II. Thus, the only real problem is the: 
determination of the potential in region MI. 

The potential in region III is not affected by the tip and 
is, therefore, the same potential as would be obtained if the 
fins were infinitely Jong. With the coordinate axes chosen 
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as shown in figure 28 (c) the point (z,y,z) is restricted to | be used to find the potential in the part of the region of 
values of y which are negative while the values of z are | interaction which is not affected by the tip. 

restricted to positive values. Note that for a tail with The velocity component normal to the fin in the z=0 
finitely long fins, the potential at a point (2,y,2) in the | plane is given by 
region of interaction is independent of the disturbances $3(%,y,0°)=—py 

produced at points located so that their projection on the 

yz-plane does not lie in the second quadrant. The general | and the velocity component normal to the fin in the y=0 
method previously derived for finding the potential resulting | Plane is given by l 

from two plane surfaces intersecting at right angles can thus i $y(x%,07,2) =pz 
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FIGURE 21.—A cross section of the surfaces of discontinuity which can be used to evaluate the potential function œz for 
two planes intersecting at 45°. 


FIGURE 22.—The addition of $, and q, to obtain the potential function œ for two surfaces intersecting at 45°. 


REPORT 11483—NATIONAL ADVISORY COMMITTEE FOR ABRONAUTICS 


zZ 


Zz 


Da 


/ 
Fiaurp 23.—A cross section of the boundary conditions of the functions $, and q for two disturbing surfaces intersecting at 60°. 


FIGURE 24.—The reflection of ¢; through one surface. 
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FIGURE 25.—A method of defining ¢; so as to eliminate discontinuities in the potential function. 
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FIGURE 26.—A method of defining q» go as to eliminate discontinuities in the potential function. 


FIGURE 27.—The addition of ¢; and $. 
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(a) One fin. (b) Two fins. Z_1,“--Mach line from 
(0) Four fins. (d) Six fins. root section 
(e) Eight fins. 
FIGURE 28.— Types of tails treated. À 
— — — Mach lines 
Figure 30 illustrates this type of normal-velocity distribution, x : 
and figure 31 illustrates a cross section of the surfaces of 7 (a) One fin. (b) Two fins, 
normal-velocity discontinuity, which previous results show (c) Four fins. (d) Six fins. 
can be used to obtain the potential in the part of the region (e) Eight fins. | 
of interaction which is not affected by the tips. FiGura 29.—Regions of similar disturbances for tails consisting of 
rectangular fins. 


Note that in figure 32 the discontinuity in the normal 
velocity across the z=0 plane is the same type as the dis- | (50) and is given as follows: 
continuity in the y=0 plane. Thus, if the potential resulting T 
from this type of discontinuity (see fig. 32) is known, then | For 03055, 
the potential resulting from any combination of discontinu- p 
ities of this type can be found. The potential for this type dlz,p,0)=— 3 [2 
of distribution is denoted by $4. By the use of cylindrical 
coordinates, as shown in figure 32, the potential at the point P?pY(1—2 cos? 6)ln Be A 
(x,p,0) can be expressed in terms of o by t+ yt — pp 


1 


ee = Bp cos 0 _ 
pleo) = holz, —4e( 2,005) (64) ANPR ae yai— B* p? sin? 0 
Equation (64) follows from figures 31 and 32. 26%? cos 0 sin 9 tan”! ra) 654 
The potential function ¢ was evaluated by use of equation P EE 4 x?— B* p? ee) 
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Por = . The potential in region IL of the fin surface is a special 
or5S05Sr, Ml 
case (03 of the preceding equation. Thus, the potential 
dot p 0) =— as E a?— ppt in region IL is given (in Cartesian coordinates) by 
; 2pz . _, BZ Bz 
29(1—2 " ln ——$e_—_ 07, = (e sin E ln ——=== } (67 
B*p ( cos ) TO ma e d(x, z) T, T Bz CN pz? ( ) 
E Bp cos O From equation (67) the pressure-difference coefficient 1s 
1 
228 p cos 0 sin Per aint ee) found to be 
g cot @ ph ae Vi dit = (68) 
28? p? sin 6 cos 6 tan! an) (65b) . 


From equations (64) and (65), the potential function in the 
region of interaction, which is not affected by the tips, is 


given by the following equation for Z ESO 


Bp 
42,0, =—La4 28*p*(1—2 cos? NR ay AY 


2xBp | cos 6 sin”? hl) 
~ \WVz?— pp? sin? 0 


i da Bp sin 0 
a ra) |* 


287p* cos 6 sin 6 (tam EOL 
+(x? — p?p? 


no! an a | | (66) 


FIGURE 31.—A cross section of the velocity discontinuity distribution 
used to find the potential in part of the region of interaction for 
a tail of four fins, 


FIGURE 30.—A cross section of the normal-velocity distribution on two 
plane surfaces representing the region of interaction for a rolling FIGURE 32.—A cross section of the velocity discontinuity distribution 
tail with four fins. associated witb the function ¢ for a tail with four fins. 


868 


As previously stated, the potential in region IV (see fig. 29) 
is a combination of the potentials in regions I, O, and IMI. 
Assume that the fins are infinitely long. In this case only 
two regions (II and III) exist since regions I and IV are 
affected by the tip. The effect of the tip can be taken into 
account by adding a potential which has zero normal velocity 
on the fin and the negative of the pressure of the infinite fin 
in the plane of the fin outboard of the tip. The value of 
such a potential on the fin is given by the difference between 
the potential of region I and the potential of region I. (This 
potential is only the effect of the tip on a semi-infinite rolling 
wing.) Thus, the potential in region 1V is the potential in 
region LU plus the difference between the potential of region 
I and the potential of region IL. Mathematically, the po- 
tential in region IV is given by 


x B UN 2) 
é(2,07,2) =P | cos”? e+e sin”? e (3—2, ne 


—— 2 — 


b z b Bz 
G=) 5) ere 
69 


From equation (69) the pressure-difference coefficient is 
given by 


TAIL CONSISTING OF SIX FINS 


The pressure and potentials on the surface of the tail con- 
sisting of six fins can be obtained in a manner similar to that 
used for the tail consisting of four fins. The pressure and 
potentials in regions I and IT (fig. 29 (d)) are the same as the 
pressure and potentials in the corresponding regions for tails 
consisting of one, two, or four fins. Regions IL and IV are 
affected by the interaction between adjacent fins. The po- 
tential in region 1V is made up of a combination of the po- 
tential in regions I, II, and 111; therefore, the main problem, 
as for the case of four fins, is the determination of the poten- 
tial in region ITI. 

The potential in region ILI is the same as the potential for 
a tail consisting of six infinitely long fins. The induced ve- 
locities normal to two of the planes of the fins are illustrated 
in figure 33. For two plane surfaces parallel to the stream 


direction and intersecting at an angle of 60°, the potential . 


in region UT can be obtained by a distribution of discontinui- 
ties in velocity as illustrated in figure 34. Note that the 
potential in region 111 can be made up of a combination of 
the potentials from a velocity discontinuity as shown in 
figure 35. The potential from this type of discontinuity is 
denoted by $. By use of cylindrical coordinates as shown 
in figure 35, the potential at the point (x,p,0) can be expressed 
in terms of dy by ; 


p(x,p,0)=40 (=,p,0-2) t (2,0,0-Z) ae (2,0047) (71) 
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Equation (71) follows from figures 34 and 35. 
The potential function ¿y was evaluated by use of equation 
(50) and is given by the following equation for 0<0<r: 


glep) PLE (x—Bo sin 6) (72) 


From equations (71) and (72) the potential function in the 
region of interaction, which is not affected by the tips, is 
T 5r, 


given by the following equation for 2 s ETE 


az,o) 22 [V3 (1—2 cos? 6)+2 cos 0 sin 8] (73) 


The potential in region IIL of the fin surface is a special 
case (9=3) of equation (73). The potential in this region is 
given (in Cartesian coordinates) by 


62,07 2) =P (74) 


z 


FIGURE 33.—A cross section of the velocity distribution normal to the 
planes of two fins of a rolling tail consisting of six fins. 
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FIGURE 34.—A. cross section of the velocity discontinuity distribution 
used to find the potential in part of the region of interaction for 
a tail of six fins. 
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From equation (74) the pressure-difference coefficient in 
region IIT is found to be zero. 

The potential in region IV is a combination of the poten- 
tials in regions 1, II, and U1; this can be shown in the same 
way as the potential in region IV of the tail consisting of 
four fins was shown to be a combination of potentials from 
other regions. Specifically, the potential in region IV for 
the tail consisting of six fins is the potential in region III 
plus the difference between the potentials of region I and of 
region IT. Mathematically, the potential in region IV is 
given by 


6 
a 2p|r2 3x2 _, p 372) 
$(x,0 Maps g C08 A 


z 22 a6) | (75) 


From equation (75) the pressure-difference coefficient is 
given by 


AG a . a E, e l 
4 rB V £ 2 2 BS 


(76) 
TAIL CONSISTING OF EIGHT FINS 


The pressure and potential on the surface of the tail 
consisting of eight fins can be found by utilizing the potential 
functions ¿y used in finding the pressure and potentials on 
the surface of the tail consisting of four and six fins. The 
pressure and potentials in regions I and II are the same as 
the pressure and potentials for the corresponding regions of 
the other tails. The potentials in regions UI, IV, V, and 
VI (see fig. 29(e)) are affected by the interaction between 
adjacent fins. Since the potentials in regions V and VI are 
combinations of the potentials in the remaining regions, the 
main problem is to find the potentials in regions III and IV. 

The potentials in regions IIT and IV are the same as the 
potential for a tail of eight infinitely long fins. The induced 
velocity normal to two of the planes of the fins is illustrated 
in figure 36. From the results for two plane surfaces inter- 
secting at an angle of 45°, the potentials in regions LIT and 
TV can be obtained by a distribution of discontinuities in 
velocity as illustrated in figure 37. The potential resulting 
from the distribution of discontinuities in velocity as illus- 
trated in figure 37 can be obtained from a distribution of 


FIGURE 35,—A cross section of the velocity discontinuity distribution 
associated with the function ¿y for a tail of six fins, 
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discontinuities as possessed by the potential function ¢y 
used in connection with the four-finned tail. The potential 
function ¿y used in connection with the four-finned tail 
was evaluated only in the region affected by the root 
sections of the fins. For the case of the eight-finned 
tail, interaction occurs between adjacent fins in regions 
which are not affected by the root sections of the fins. 
The potential function ¢ of the four-finned tail in the region 
not affected by the root section must be known. In this 
region, the potential functions ¢ for the four- and six-finned 
tails are the same. 

From figures 32 and 37 the potential function in the region 
of fin interaction which is affected by the root sections of 


the fins is, for I<0sto 


d(x, p, = — polz, 0,9) + do (2.0.0 m 


do (2.0.0) +4 (2,0,045) (77) 


FIGURE 36.—A cross section of the normal velocity induced on the 
planes of two fins of a rolling tail consisting of eight fins. 


FIGURE 37.—A cross section of the velocity discontinuities which can 
be used to obtain the potential in the region of interaction between 
adjacent fins for a tail consisting of eight fins. 


870 ‘ REPORT 1143—NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


where ¢ is given by equations (65). 
The ae in the part of the region of interaction which is affected by the root section is (from eqs. (65) and (77)), 


for 25055) 
Bp cos (9-3) P 
re) ES —z cos O sin” eat? COs (0-4) sin! 4 ee - Bp = 0 
TB x 2 — B p ?— Bp? sin? (0-7) yT — B*p cos? 0 


re eoe ae a GP 
Ae (0-7) Es TEP | (78) 


Bp 3 
aroso mnt Ja e epe 


The potential in region IL is a special case (93) of equation (78). Setting => in equation (78) yields (in Cartesian 
coordinates) 


= _ 2pz ==" _ —] 2 ws wv 
p(1,07,2)= =B (sin E 2/2 sin eat Pe tan ea (79) 
From equation (79) the pressure-difference coefficient is found to be given by 
, _ 8p2 f. BZ „a Bz 
AC ey sin”! —-—y2 sin Tr pa (80) 


By inspection of figures 35 and 37 the potential function in the region of fin interaction, which is not affected by the 


root sections of the fins, is, for 5 <0 sz 
o(p) =t (20,0)-+40 (2,0,0—F)—60 (2,0.05)—60 (2100) (81) 


where ¢ is given by equation (72). Substituting equation (72) into equation (81) yields (remember that | de 18 Zero 
upstream of the Mach cone from the y-axis), for gts” l 


(x, p 0) =P =| (1-7) sin sin oe oe [or cos 6 sin 922 (1—2 sin? 0)— 
pct ados 
PP (cos0+sin0) 8 E cos} l , ° (82) 
By/2 E (sino—cos0)ita > 2 (in0— —cosó) B z—Bp sin 0; if x>Bp sin 0 


The potential in region IV is a special case (= 2) of equation (82). Setting => in equation (82) yields (in Cartesian 
coordinates) 


a 


$(z, 07, y [G— /2)2+82] (83) 
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From equation (83) the pressure-difference coefficient is found to be given by 


_ Ape, i 


The potential in region V is the potential in region IV plus the difference between the potential in region I and the 
potential in region IJ; thus, from equations (62a), (62b), and (83), the potential in region V is found to be given by 


b | 
#2, 0- y 22(.21D 422 sn" 5) 1, y) EAE 88) 


From equation (85) the pressure-difference coefficient is found to be given by 


AC e. sin! e). NE! | (86) 
”  *"BVLJa \ Y V 2 2'6 l 


The potential in region VI is the potential in region IO plus the difference between the potential in region I and the 
potential in region I; thus, from equations (62a), (62b), and (79), the potential in region VI is found to be given by 


b 
ae eee tan! —_*__+-x2 sin”! EN) Cual 
y 21 — p222 yxi—p8?2? t 


EA o 


From equation (87), the pressure-difference coefficient in region VI is found to be given by 


| T) 
___8p -1 $Z ee _ o(3-=) b JE 88 
somi: A O C = 


DISCUSSION OF RESULTS FOR ROLLING TAILS 


- na 2pz| _ abe. ss 
olz, 0 azl % COS”! -—2 2 sin 


Illustrative plots of the chordwise and spanwise pressure distributions across one fin for tails with various numbers of 
fins are shown in figure 38. Figure 39 shows illustrative plots of the spanwise loadings on one fin for tails made up of 
various numbers of fins. l 

The potential function $y used in finding the pressures and ‘potentials for the tails consisting of four, six, and eight fins 
could be used in finding pressures and potentials for tails consisting of any even number of fins provided that the region of 
interaction between adjacent fins does not affect the tip. The restriction on the region of interaction causes the range of 
validity to decrease as the number of fins is increased. The range of validity could be extended, however, by use of a pressure 
or potential cancellation method such as given in references 14 and 15. 

From the potential, the damping in roll per fin was calculated. Table I presents the results of these calculations. Figure 
40 presents the variation of the damping in roll per fin with 48 for tails made up of various numbers of fins. Fora given 
Mach number (8 constant), figure 40 shows the variation of the damping with aspect ratio. Figure 41 presents the variation 
of the damping in roll per fin with Mach number for tails consisting of various numbers of fins with a fin aspect ratio of 1.5. 


LANGLEY ABRONAUTICAL LABORATORY, 
NATIONAL Apvisory COMMITTEE FOR AERONAUTICS, 
Lanauey Frer, Va., October 26, 1951. 
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(a) Tail of one fin. 


FiGurE 38.—Chordwise and spanwise pressure distributions on a fin of aspect ratio 1.5 at a Mach number of y2 for tails 
consisting of rectangular fins. 
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(b) Tail of two fins. 
Freurs 38.—Continued. 
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cal Mach lines 


(c) Tail of four fins. 
FiGURB 38.—Continued. 
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(d) Tail of six fins. 
FIGURE 38.—Continued. 
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(e) Tail of eight fins. 
FIGURE 38.—Concluded. 


—  2y/b 
Fiaurb 39.—Spanwise loading on a fin of aspect ratio 1.5 at a Mach 
number of y2 for tails consisting of rectangular fins, 
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FIGURE 40.—The damping in roll per fin for tails consisting of rectangular fins. 
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FIGURE 41,—An illustrative variation of Ci, per fin with Mach number for tails consisting of rectangular fins with a fin aspect ratio of 1.5 
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TABLE I 


` DAMPIN G-IN-ROLL COEFFICIENT PER FIN 


BC,, per fin 


Valid for 
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